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Abstract 

Computing Delaunay triangulations in involves evaluating the so-called in_sphere predi¬ 
cate that determines if a point x lies inside, on or outside the sphere circumscribing d+ 1 points 
Po 7 ... ,pd. This predicate reduces to evaluating the sign of a multivariate polynomial of degree 
d + 2 in the coordinates of the points x,po,... ,pd- Despite much progress on exact geometric 
computing, the fact that the degree of the polynomial increases with d makes the evaluation 
of the sign of such a polynomial problematic except in very low dimensions. In this paper, we 
propose a new approach that is based on the witness complex, a weak form of the Delaunay 
complex introduced by Carlsson and de Silva. The witness complex Wit(L, W) is defined from 
two sets L and W in some metric space A': a finite set of points L on which the complex is 
built, and a set W of witnesses that serves as an approximation of A'. A fundamental result of 
de Silva states that Wit(A, W) = Del(L) if W = X = R d . In this paper, we give conditions on 
L that ensure that the witness complex and the Delaunay triangulation coincide when W is a 
finite set, and we introduce a new perturbation scheme to compute a perturbed set L close to 
L such that Del(L') = Wit(L , ,IT). Our perturbation algorithm is a geometric application of 
the Moser-Tardos constructive proof of the Lovasz local lemma. The only numerical operations 
we use are (squared) distance comparisons (i.e., predicates of degree 2). The time-complexity 
of the algorithm is sublinear in \W\. Interestingly, although the algorithm does not compute 
any measure of simplex quality, a lower bound on the thickness of the output simplices can be 
guaranteed. 
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1 Introduction 


The witness complex was introduced by Carlsson and de Silva (dSC04| as a weak form of the 
Delaunay complex that is suitable for finite metric spaces and is computed using only distance 
comparisons. The witness complex Wit(L, W) is defined from two sets L and W in some metric 
space X: a finite set of points L on which the complex is built, and a set W of witnesses that serves 
as an approximation of X. A fundamental result of de Silva |dS08] states that Wit(L, W) = Del(L) 
if W is the entire Euclidean space X = M d , and the result extends to spherical, hyperbolic and 
tree-like geometries. The result has also been extended to the case where W = X is a smoothly 
embedded curve or surface of M. rf [AEM07J. However, when the set W of witnesses is finite, the 
Delaunay triangulation and the witness complexes are different and it has been an open question 
to understand when the two structures are identical. In this paper, we answer this question and 
present an algorithm to compute a Delaunay triangulation using the witness complex. 

We first give conditions on L that ensure that the witness complex and the Delaunay trian¬ 
gulation coincide when W C is a finite set (Section |3|. Some of these conditions are purely 
combinatorial and easy to check. In a second part (Section [I]), we show that those conditions can 
be satisfied by slightly perturbing the input set L. Our perturbation algorithm is a geometric appli¬ 
cation of the Moser-Tardos constructive proof of the general Lovasz local lemma. Its analysis uses 
the notion of protection of a Delaunay triangulation that we have previously introduced to study 
the stability of Delaunay triangulations |BDG13| . 

Our algorithm has several interesting properties and we believe that it is a good candidate for 
implementation in higher dimensions. 


1. Low algebraic degree. The only numerical operations used by the algorithm are (squared) 
distance comparisons (i.e., predicates of degree 2). In particular, we do not use orientation or 
in-sphere predicates, whose degree depends on the dimension d and are difficult to implement 
robustly in higher dimensions. 


2. Efficiency. Our algorithm constructs the witness complex Wit(L / , W) = Del(L') of the per¬ 

turbed set L' in time sublinear in \W\. See Section [5] 

3. Simplex quality and Delaunay stability. Differently from all papers on this and related 

topics, we do not compute the volume or any measure of simplex quality. Nevertheless, 
through protection, a lower bound on the thickness of the output simplices can be guaranteed 
and the resulting Delaunay triangulation is stable with respect to small 


(see Theorem 26), 
metric or point perturbations |BDG13j . 


4. No need for coordinates. We can construct Delaunay triangulations of points that come from 
some Euclidean space but whose actual positions are unknown. We simply need to know the 
interpoint distances. 


5. A thorough analysis. Almost all papers in Computational Geometry rely on oracles to eval¬ 
uate predicates exactly and assume that the complexity of those oracles is 0(1). Our (prob¬ 
abilistic) analysis is more precise. We only use predicates of degree 2 (i.e. double precision) 
and the analysis fully covers the case of non generic data. 
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1.1 Previous work 


Millman and Snoeyink jMSIO] developed a degree-2 Voronoi diagram on a U x U grid in the 
plane. The diagram of n points can be computed using only double precision by a randomized 
incremental construction in 0{n log n log U) expected time and O(n) expected space. The diagram 
also answers nearest neighbor queries, but it doesn’t use sufficient precision to determine a Delaunay 
triangulation. 

Our work has some similarity with the e-geometry introduced by Salesin et al. |SSG89j . For a 
given geometric problem, the basic idea is to compute the exact solution for a perturbed version of 
the input and to bound the size of this implicit perturbation. The authors applied this paradigm 
to answering some geometric queries such as deciding whether a point lies inside a polygon. They 
do not consider constructing geometric structures like convex hulls, Voronoi diagrams or Delaunay 
triangulations. Also, the required precision is computed on-line and is not strictly limited in advance 
as we do here. 

Further developments have been undertaken under the name of controlled perturbation jHallOj . 
The purpose is again to actually perturb the input, thereby reducing the required precision of 
the underlying arithmetic and avoiding explicit treatment of degenerate cases. A specific scheme 
for Delaunay triangulations in arbitrary dimensions has been proposed by Funke et al. |FKMS05] , 
Their algorithm relies on a careful analysis of the usual predicates of degree d + 2 and is much more 
demanding than ours. 

1.2 Notation 

In the paper, W and L denote sets of points in T d where T d denotes the standard flat torus M d /Z d . 
L is finite but we will only assume that W is closed in T d . The points of W are called the witnesses 
and the points of L are called the landmarks. To keep the exposition simpler, boundary issues will 
be handled in the full version of the paper (see Section [6]). 

We say that W C T d is an e-sample if for any x G T d there is a w G W with ||u; — ,t|| < e. 
If L C T d is a finite set, then it is a A'-sample for T d for all X' greater or equal to some finite 
A > 0. The parameter A is called the sampling radius of L. We further say that L is (A ,p)-net if 
\\p ~ Q || > fi A for all p,q G L. We call p the sparsity ratio of L. Note that, for any two points p, q 
of a ( X,p)-net , we have pX < \\p — q\\ < 2A, which implies p < 2. 

A simplex a C L is a finite set. We always assume that L contains a non-degenerate d-simplex, 
and we demand that A < 1/4. This bound on A avoids the topological complications associated 
with the periodic boundary conditions. 

The volume of Euclidean d-ball with unit radius will be denoted by Ud- 


2 Delaunay and witness complexes 

Definition 1 (Delaunay center and Delaunay complex) A Delaunay center for a simplex a C 
L is a point x G T d that satisfies 

||p — ®|| < \\q — a:||, Vp G a and Vg G L. 

The Delaunay complex Del(L) of L is the complex consisting of all simplexes a C L that have a 
Delaunay center. 
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Figure 1: Blue points are the witnesses W and the black points are the landmarks L. The simplex 
a is witnessed by w. 


Note that x is at equal distance from all the vertices of a. A Delaunay simplex is top dimensional 
if is not the proper face of any Delaunay simplex. The affine hull of a top dimensional simplex has 
dimension d. If u is top dimensional, the Delaunay center is the circumcenter of a which we denote 
c a . We write R a for the circumradius of a. 

Delaunay |Del34] showed that if the point set L is generic, i.e., if no empty sphere contains d+ 2 
points on its boundary, then Del(L) is a triangulation of T d (see the discussion in Section [ 3 ]), and 
any perturbation L' of a finite set L is generic with probability 1. We refer to this as Delaunay’s 
theorem. 

We introduce now the witness complex that can be considered as a weak variant of the Delaunay 
complex. 

Definition 2 (Witness and witness complex) Let a be a simplex with vertices in L C T d , and 
let w be a point of W C T d . We say that w is a witness of a if 

||u; — p\\ < ||w — g||, Vp E a andVq£L\(T. 

The witness complex Wit(L, W ) is the complex consisting of all simplexes a such that for any simplex 
t C a, t has a witness in W. 

Observe that the only predicates involved in the construction of Wit(L, W) are (squared) dis¬ 
tance comparisons, i.e. polynomials of degree 2 in the coordinates of the points. This is to be 
compared with the predicate that decides whether a point lies inside, on or outside the sphere 
circumscribing a d-simplex which is a polynomial of degree d + 2. 

3 Identity of witness and Delaunay complexes 

In this section, we make the connection between Delaunay and witness complexes more precise. We 
start with de Silva’s result |dS08j : 

Theorem 3 Wit(L, W) C Del(L), and ifW = T d then Wit(L, W) = Del(L). 
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If L is generic, we know that Del(L) is embedded in T d by Delaunay’s theorem. It therefore 
follows from Theorem [3] that the same is true for Wit(L,VF). In particular, the dimension of 
Wit(L,W) is at most d. 


3.1 Identity from protection 

When W is not the entire space T d but a finite set of points, the equality between Del(L) and 
Wit(L,W) no longer holds. However, by requiring that the d-simplices of Del(L) be 5-protected, 
a property introduced in [BDG13| . we are able to recover the inclusion Del(L) C Wit(L, W), and 
establish the equality between the Delaunay complex and the witness complex with a discrete set 
of witnesses. 


Definition 4 (d-protection) We say that a simplex a C L is 5-protected at x € T d if 

\\x — g|| > ||x — p\\ + 5, Vp 6 o and Vq £ L \ a. 

We say that Del(L) is 5-protected when each Delaunay d-simplex of Del(L) has a d-protected 
Delaunay center. In this sense, d-protection is in fact a property of the point set and we also say 
that L is d-protected. If Del(L) is d-protected for some unspecified d > 0, we say that L is protected 
(equivalently L is generic). We always assume d < A since it is impossible to have a larger d if L is 
a A-sample. The following lemma is proved in |BDG014| . For a subcomplex K 6 Del(L), we define 
star(A'; Del(L)) as the subcomplex consisting of all the simplices that have at least one vertex in 
K (note that this definition departs from common usage). We will use star 2 (p) as a shorthand for 
star(star(p; Del(T)); Del(L)). 

Lemma 5 (Inheritance of protection) Let L be a ( X,fi)-net and suppose p £ L. If every d- 
simplex in star 2 (p) is 5-protected, then all simplices in star(p; Del(L)) are at least 5'-protected where 
5' = Specifically, any Delaunay k-simplex r is 5’-protected at a point z T which is the barycenter 
of the circumcenters of a subset of d — k + 1 d-cofaces of r. 

The following lemma is an easy consequence of the previous one. 


Lemma 6 (Identity from protection) Let L be a ( X,p,)-net with p E L. If all the d-simplices 
in star 2 (p) are 5-protected and W is an e-sample for T d with 5 > then star(p; Wit(L, W)) = 
star(p; Del(L)). 


Proof By Theorem [3j we have 

star(p; Wit(L, W)) C star(p; Del(L)). 

We now prove the other inclusion. Let a be a simplex in star(p; Del(L)) and set 

g = ^>2e. 

4 d 

Let r C a be a face of a. Then, by Lemma [5] t is d'-protected at a point z r such that 

1. \\z T - Pi\\ = \\z T - pj\\ = r V Pi , Pj £ t 

2. \\z T — pi\\ > r + 5' Mpi G L\t 
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For any x E B(z T , <5'/2), any pi E r and any pi £ L\t, we have 


Ik -Pill < Ikr - Pill + I \z T - x\\ < r + — 

and 

„ .. < 5 ' 

Ik - Pill > Ikr — Pill - Ik — z t\\ > r + o - — = r + —. 

Hence, x is a witness of r. Since e < 5'/ 2, there must be a point re E W in B(z T ,5'/ 2) which 
witnesses r. Since this is true for all faces r C cr, we have cr E star(p; Wit(L, IT)). □ 


We end this subsection with a result proved in |BDG13i Lemma 3.13] that will be useful in 
Section [5] For any vertex p of a simplex a, the face oppposite p is the face determined by the other 
vertices of a, and is denoted by a p . The altitude of p in a is the distance D(p, a) = d(p, aff(cq,)) from 
p to the affine hull of a p . The altitude D(a) of cr is the minimum over all vertices p of cr of D(p, cr). 
A poorly-shaped simplex can be characterized by the existence of a relatively small altitude. The 
thickness of a j-simplex a is the dimensionless quantity ©(cr) that evaluates to 1 if j = 0 and to 
jA(a) °therwise, where A(cr) denotes the diameter of cr, i.e. the length of its longest edge. 


Lemma 7 (Thickness from protection) Suppose a E Del(L) is a d-simplex with circumradius 
less than A and shortest edge length greater than or equal to JxX. If every (d— 1 )-face of a is also a 
face of a 5-protected d-simplex different from a, then the thickness of a satisfies 


Q(a) > 


5 (fl + d) 

8 d 


In particular, suppose p E L, where L is a (A ,p)-net, and every d-simplex in star 2 (p) is 5- 
protected, then every d-simplex in star(p) is (§§) -thick. 


3.2 A combinatorial criterion for identity 


The previous result will be useful in our analysis but does not help to compute Del(L) from 
Wit (L,W) since the ^-protection assumption requires knowledge of Del(L). A more useful re¬ 
sult in this context will be given in Lemma 10 Before stating the lemma, we need to introduce 


some terminology and, in particular, the notion of good links. 

A complex K is a k-pseudo-manifold if it is a pure ^-complex and every (k — l)-sinrplex is the 
face of exactly two L’-simplices. 


Definition 8 (Good links) Let K be a complex with vertex set L C T d . We say p E L has a good 
link i/link(p;ii) is a (d — l)-pseudo-manifold. If every p E L has a good link, we say K has good 
links. 

For our purposes, a simplicial complex K is a triangulation of T d if it is a d-manifold embedded in 
T d . We observe that a triangulation has good links. 


Lemma 9 (Pseudomanifold criterion) If K is a triangulation of and J C K has the same 
vertex set, then J = K if and only if J has good links. 
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Figure 2: Pinched simplicial torus is an example of 2-pseudo-manifold. 


Proof Let L be the common vertex set of K and J, and let p be a point in L. Since K is 
a triangulation, link(p; K ) is a simplicial (d — l)-sphere, which is a manifold. A (d — 1 ^pseudo- 
manifold Cj cannot be properly contained in a connected (d — l)-manifold simplicial complex Ck, 
because a (d — l)-simplex that does not belong to Cj cannot share a (d — 2)-face with a simplex 
of Cj\ Every (d — 2)-simplex of Cj is already the face of two (d — l)-simplices of Cj. Since any 
two (d — l)-simplices in Ck can be connected by a path confined to the relative interiors of (d — 1) 
and (d — 2)-simplices, it follows that a (d — l)-simplex r E Ck must also belong to Cj. For an 
arbitrary simplex a E Ck, we have that a is a face of a (d — l)-simplex r E Ck, because Ck is a 
pure (d — l)-complex (it is a manifold). It follows that r E Cj and therefore a E Cj. 

Therefore, since link(p; J) C link(p; K), if J has good links we must have link(p; J) = link(p; K ). 
It follows that star(p; J) = star(p; K ) for every p E L, and therefore J = K. If J does not have 
good links it is clearly not equal to K , which does. □ 

We can now state the lemma that is at the heart of our algorithm. It follows from Theorem [3j 
Lemma [9j and Delaunay’s theorem: 

Lemma 10 (Identity from good links) If L is generic and the vertices o/Wit(L, W) have good 
links, then Wit(L, IT) = Del(L). 


4 Turning witness complexes into Delaunay complexes 


Let, as before, L be a finite set of landmarks and W a finite set of witnesses. In this section, 


we intend to use Lemma 10 to construct Del(d/), where L' is close to L, using only comparisons 


of (squared) distances. The idea is to first construct the witness complex Wit(L,lT) which is a 
subcomplex of Del(L) (Theorem [3]) that can be computed using only distance comparisons. We 
then check whether Wit(L, W) = Del(L) using the pseudomanifold criterion (Lemma [9]). While 
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there is a vertex p of Wit(L, W) that has a bad link (i.e. a link that is not a pseudomanifold), we 
perturb p' and the set of vertices I(p'), to be exactly defined in Section 4.3 (See Eq. [TJ that are 
responsible for the bad link L(p') = link(j/, Wit(L', W), and recompute the witness complex for 
the perturbed points. We write L' for the set of perturbed points at some stage of the algorithm. 
Each point p' is randomly and independently taken from the so-called picking ball B(p,p). Upon 
termination, we have Wit(L 7 , W) = Del(L'). The parameter p, the radius of the picking balls, must 
satisfy Eq. Q to be presented later. The steps are described in more detail in Algorithm [l] The 
analysis of the algorithm relies on the Moser-Tardos constructive proof of Lovasz local lemma. 


Algorithm 1 Delaunay triangulation from witness complex 
Input: L, W, p , A, p 

while a vertex p' of Wit (I/, W) has a bad link do 
perturb p' and the points in I(p') (defined in Eq. [I]) 
update Wit(L', W) 

end while 

Output: Wit(L', W) = Del(L') 


4.1 Lovasz local lemma 

The celebrated Lovasz local lemma is a powerful tool to prove the existence of combinatorial ob¬ 
jects |AS08] , Let A be a finite collection of “bad” events in some probability space. The lemma 
shows that the probability that none of these events occur is positive provided that the individual 
events occur with a bounded probability and there is limited dependence among them. 


Lemma 11 (Lovasz local lemma) Let A = {A\, ..., Ajy} be a finite set of events in some prob¬ 
ability space. Suppose that each event Ai is independent of all but at most T of the other events Aj, 
and that Pr [Aj] < vo for all 1 < i < N. If 

m -ArTT) 

(e is the base of the natural logarithm), then 


Pr 


‘ N 


_i=1 


> o. 


Assume that the events depend on a finite set of mutually independent variables in a probability 
space. Moser and Tardos |MT10| gave a constructive proof of Lovasz lemma leading to a simple 
and natural algorithm that checks whether some event A E A is violated and randomly picks new 
values for the random variables on which A depends. We call this a resampling of the event A. 
Moser and Tardos proved that this simple algorithm quickly terminates, providing an assignment of 
the random variables that avoids all of the events in A. The expected total number of resampling 
steps is at most N/T. 


4.2 Three geometric lemmas 


We will be using the following geometric lemmas. 


We will need the following technical result to prove Lemma 15 
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Figure 3: For Lemma 


14 


Lemma 12 If L is a \-sample of T d , the circumradms R& of any simplex a in Wit (L,W) is less 
than X. The witnesses of a are at a distance less than A from the closest vertex of a and at a 
distance less than 3A from any vertex of a. 

Proof By Theorem [ 3 ] a is a simplex of Del(L) and, since L is a A-sample of T d , the circumradius 
R a is less than A. Let re be a witness of 0 and q the vertex of o closest to w. Since L is a A-sample, 
||rc — q || < A. The diameter of a is thus at most 2 R a < 2A, and ||u; — p|| < 3A for any vertex p E 0 . 

□ 


Lemma 13 If L is a (A ,p)-net ofT d and 4 p < p, then L' is an (A ',p')-net, where 

X' = A(1 + p) and p = -- 

1 + p 3 


Lemma 14 (See [BDG14L Lemma 5.2]) Let S(c,R ) be a hypersphere o/M rf of radius R centered 
at c and T$ the spherical shell = B(c,R + 5) \ B(c, R). Let, in addition, B p denote any d-ball of 
radius p < R. We have 

vol d (T s n B p ) < U d - 1 (^p^j 6, 

where U d _\ denotes the volume of a Euclidean (d — 1 )-ball with unit radius. 

4.3 Correctness of the algorithm 

We first bound the density radius A 7 and the sparsity ratio p' of any perturbed point set L'. We 
will write p = pX and p = pX, and assume p < p/ 4. 

We refer to the terminology of the Lovasz local lemma. Our variables are the points of L' which 
are randomly and independently taken from the picking balls B(p, p), p E L. 
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The events are associated to points of L ', the vertices of Wit (L',W). We say that an event 
happens at p' £ L when the link L(p') of p' in Wit(L / , W) is not good, i.e., is not a pseudomanifold. 
We know from Lemma [h] that if p’ is a vertex of Wit(L / , W) and L(p') is not good, then there must 
exist a d-simplex in star 2 ip') that is not d-protected for 5 = 8 de/p,. We will denote by 


• hip') : the set of points of L' that can be in star 2 (p') = star(star(j/; Del(L'); Del(L')) 

• I 2 ip') : the set of points of L' that can violate the d-protected zone Zs(a') = B(c a , R a + d) \ 
B(c a ,R a ) for some d-simplex a' in star 2 (p') 

• I{p') : Since computing I\ (p r ) and exactly is difficult, we will rather compute the set 


I ip') = { q' ■ q £ L n B(p, R) where R = 5A + | 


(1) 


Using triangle inequality, Lemma 13, and the facts that p < p /4 and 5 < A, we can show that 

I{p) 2 hip) U/ 2 (pO> 

where Ii{p') is the set of points in L that correspond to points in hip') C L'. 

Sip') : the set of d-simplices with vertices in hip') that can belong to star 2 (y/) 


The probability w\ ip') that L(pf) is not good is at most the probability W 2 ip ') that one of the 
simplices of Sip'), say a', has its d-protecting zone Zsicr') violated by some point of l!. Write 
w^iq',a') for the probability that q' belongs to the d-protection zone of the d-sinrplex a'. We have 


-uJlip') <-w 2 ip') < w 3iq',h) ( 2 ) 

q'£h(p') o-'eS(p') 

The following lemmas provide upper bounds on \I{p')\, |<S'(p , )|, T and w^iq', a'). 


Lemma 15 |/(p')| < / = ( h ) and, |S'(p , )| < K= 


jd+1 

WIL 


Proof Since any two points in L are at least pX apart, a volume argument then gives the announced 
bound. Referring to 0 we have 


\Hp')\ < 


'5A+ ^ 

jl A 

, ~T 


10 


< 4 + - < hr 


We thus have 


W)\ < 


< 


p 


jd +1 


14 


P 


Lemma 16 An event is independent of all but at most T 


d + 1 J (d + 1)! 


□ 


p J 


other events. 
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Proof Two events p' and q' are independent if I(p') D I(q') = 0. Referring to R defined in Q, 
this implies p' and q' are independent if \\p — g|| < 2R = 10A + 3p\. Hence an upper bound on the 
number of events that may not be independent from the event at p' is given by a volume argument 
as in the proof of Lemma 15 Here we count points since the events are associated to points: 


10 A + 3pX + 


fi A 
AT 


< 


„ 20 
7+ — 
P 


<i” 

P 


def p 


□ 


Lemma 17 m 3 (q' , a') < 2n d 1 ^. 


Proof The volume of Euclidean d-ball with unit radius will be denoted by Ud. The probability 
w 3 (q r , a') is the ratio between the volume of the portion of the picking region that is in the spherical 
shell B(c a ',R a i + 5) \ B(c a i,R a i) and the volume of the picking region. Hence, using |BDG141 
Lem. 5.2] and the crude bound < 2 d , we obtain 


w 3 (q',a') < 


Ud-ij^pY^l 

u d p d 


< 2n d ~ 1 


5_ 

P 


□ 


Using Eq. ([2]) , and Lemmas 15 and 17, we conclude that 

^l(pO < 27T d_1 


IK 

P 


An event depends on at most T other events. Hence, to apply the Lovasz Local Lemma HU i1; 
remains to ensure that w\ (p') < In addition, we also need that 5 > yf to be able to apply 

Lemma [6] We thefore need to satisfy the following inequality: 


8 de 


~rr < S < Jp where J 1 = f 2e7r d l IK(T + 1) = I ^ ) 
P \P) 


d -1 ; 


2 A 


0{d?) 


(3) 


Observe that I, K, T and J depend only on p and d. We conclude that the conditions of the Lovasz 
local lemma hold if the parameter p satisfies 


n 24 de 

- > p > - 

4 - H - fid 


(4) 


Hence, if e is sufficiently small, we can fix p so that Eq. Q holds. The algorithm is then 
guaranteed to terminate. By Lemma 10 the output is De^Z/). 

It follows from Moser-Tardos theorem that the expected number of times the "while-loop" in 
Algorithm jlj runs is O an< I since \I(p')\ < I, we get that the number of point perturbations 

performed by Algor ithmjlj is O = 0(\L\) on expectation, where the constant in the O depends 

only on d and not on the sampling parameters A or p. We sum up the results of this section in 


Theorem 18 Under Eq. &■ Algorithm [7] terminates and outputs the Delaunay triangulation of 
some set L' whose distance to L is at most p. The number of point perturbations performed by the 
algorithm is O (|L|). 
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5 Sublinear algorithm 


When the set ll is generic, K = Wit(L 7 , W) is embedded in T d and is therefore d-dimensional. It is 
well known that the d-skeleton of Wit(L 7 , W) can be computed in time 0((| W| + |7f|) log \L'\) using 
only distance comparisons |BM14| . Although easy and general, this construction is not efficient 
when W is large. 

In this section, we show how to implement an algorithm with execution time sublinear in \W\. 
We will assume that the points of W are located at the centers of a grid, which is no real loss of 
generality. The idea is to restrict our attention to a subset of W, namely the set of full-leaf-points 
introduced in Section 5.2 These are points that may be close to the circumcentre of some d-simplex. 
A crucial observation is that if a d-simplex has a bounded thickness, then we can efficiently compute 
a bound on the number of its full-leaf-points. This observation will also allow us to guarantee some 
protection (and therefore thickness) on the output simplices, as stated in Theorem 26 below. 

The approach is based on the relaxed Delaunay complex, which is related to the witness complex, 
and was also introduced by de Silva |dS08| . We first introduce this, and the structural observations 
on which the algorithm is based. 


5.1 The relaxed Delaunay complex 

The basic idea used to get an algorithm sublinear in \W\ is to choose witnesses for d-simplices 
that are close to being circumcentres for these simplices. With this approach, we can in fact avoid 
looking for witnesses of the lower dimensional simplices. The complex that we will be computing is 
a subcomplex of a relaxed Delaunay complex: 


Definition 19 (Relaxed Delaunay complex) Let a C L' be a simplex. An a-center of a is any 
point x E T d such that 

\\x — p\\ < \\x — q\\ + a, Vp, q E cr. 

We say that x is an a-Delaunay center of a if 

||® — P || < ||® — <711 + a, Vp E a and\/q E L'. 

The set of simplices that have an a-Delaunay centre in W is a simplicial complex, called the a- 
relaxed Delaunay complex, and is denoted Del“(L 7 ,IT). 

We say w E W is an a-witness for a C L' if 

||w — p|| < ||u> — g|| + a, Vp E a and \/q E L' \ a. 

We observe that roElkis an a-Delaunay center if and only if it is an a-center and also an a-witness. 


Lemma 20 The distance between an a-Delaunay center for a E Del“(L 7 ,IT) and the farthest vertex 
in a is less than A ' + a. In particidar, < 2A 7 + 2a. 

I/rE Del(L 7 ) and c is a Delaunay center of t, then any point in B(c, r) is a 2r-Delaunay centre 
for t. Thus Del(L 7 ) C Del 2e (L 7 , W). 


If, for some 6 > 0, all the d-simplices in Del“(L 7 ,lT) have a d-protected circumcentre, then we 
have that Del“(L 7 ,lT) C Del(L 7 ), and with a > 2e, it follows (Lemma 201 that Del Q (L 7 ,lT) = 
Del(L 7 ), and Del(L 7 ) is itself d-protected and has good links. 
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Reviewing the analysis of the Moser-Tardos algorithm of Section 4.3 we observe that the exact 
same estimate of W 2 (//) that serves as an upper bound on the probability that one of the simplices 
in star 2 (p 7 , Del(L 7 )) is not 5-protected at its circumcenter, also serves as an upper bound on the 
probability that one of the simplices in star 2 (p 7 , Del“(L 7 , T d )) is not 5-protected at its circumcenter, 
provided that 4a + <5 < A 7 (using the diameter bound of 2 A 7 + 2 a from Lemma |20|) , which we will 
assume from now on. We can therefore modify Algorithm [I] by replacing Wit(L 7 , IT) by Del“(L 7 , IT). 

We now describe how to improve this algorithm to make it efficient. For our purposes it will be 
sufficient to set a = 2e. In order to obtain an algorithm sublinear in |1T|, we will not compute the 
full Del 2 ~(L 7 , W) but only a subcomplex we call Delo £ (L 7 , IT). The exact definition of Delo £ (L 7 , IT) 
will be given in Section |5.3| but the idea is to only consider d-simplices that show the properties 
of being ©o-thick for some parameter ©o to be defined later (Eq. [5] and [6]). This will allow us to 
restrict our attention to points of IT that lie near the circumcentre. As explained in Section 5.3 
this is done without explicitly computing thickness or circumcentres. 

As will be shown in Section 5.4 (Lemma 24), the modification of Algorithm [I] that computes 
Delg £ (L',W) instead of Wit(L 7 ,IT) will terminate and output a complex Delg £ (.L 7 , IT) with good 
links. However, this is not sufficient to guarantee that the output is correct, i.e., that Del 2£ (L 7 , IT) = 
Del(L'). In order to obtain this guarantee, we insert an extra procedure checkQ, which, without 
affecting the termination guarantee, will ensure that the simplices of Delo £ (L 7 , IT) have e)*-protected 
circumcenters for a positive 5*. It follows that Del q £ (L 7 ,IT) C Del(L 7 ) and, by Lemma [ 9 ] that 
Delo £ (L / , W) = Del(L 7 ). The modification of Algorithm[l]that we will use is outlined in Algorithm]^] 


Algorithm 2 Protected Delaunay triangulation from Delg £ (L',W) 


input: L, IT, p, s, A, p 
L' <- L 


compute: Delg E (L / ,VF) 

while a vertex p' of Delg^L 7 , IT) has a bad link or check(p 7 ) = FALSE do 
perturb p' and the points in I(p') ( I(p') is defined in Section [4.3| ) 

update Delo £ (L', W) 

end while 

output: DelQ £ (L / ,lT) = Del(L') with <5* protection (See Lemma 25) 


We describe the details of computing Delg £ (L / , IT) and of the check() procedure in the following 
subsections. 

5.2 Computing relaxed Delaunay centers 

We observe that the a-Delaunay centers of a d-simplex cr are close to the circumcenter of cr, provided 
that cr has a bounded thickness: 

Lemma 21 (Clustered a-Delaunay centers) Assume that L' is a (A 7 , ft')-sample. Let a be a 
non degenerate d-simplex, and let x be an a-center for a at distance at most C A 7 from the vertices 
of a, for some constant C > 0. Then x is at distance at most from the circumcenter c a of a. 

In particular, if x is an a-Delaunay center for a, then 

2 a 

\\Ca ~ X < ■ 

0(o-)p 7 
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Proof Using |IBDG131 , Lemma 4.3] and A(<r) > p'X', we get 

CX'a Ca 


\w, T - cJ < 


< 


0(cr)A(cr) @(a)p' 


The second assertion follows from this together with Lemma 20, and the assumption that a < X'. 

□ 


It follows from Lemma 21 that a-Delaunay centers are close to all the bisecting hyperplanes 

Hpq = {x G I \\x - p\\ = ||x - q\\, p,q G a}. 

The next simple lemma asserts a kind of qualitative converse: 

Lemma 22 Let a be a d-simplex and H pq be the bisecting hyperplane of p and q. A point x that 
satisfies d(x , H pq ) < a, for any p, q G a is a 2a-center of a. 

Let a be a d-simplex of Del“(L / ,VU) and let LI be the smallest box with edges parallel to the 
coordinate axes that contains a. Then the edges of Ll have length at most 2A 7 + 2a (Lemma 20). 


Any a-Delaunay center for a is at a distance at most A' + a from Ll. Therefore all the a-Delaunay 
centers for a lie in an axis-aligned hypercube LI with the same center as U and with side length at 
most 4A 7 + 4a < 5X'. Observe that the diameter (diagonal) z of LI is at most 5A '\fd. 

Our strategy is to first compute the a-centers of a that belong to Ll D W and then to determine 
which ones are a-witnesses for a. Deciding if an a-center is an a-witness for a can be done in 


constant time since L' is a (A 7 ,/r 7 )-net and A CT < 2A 7 + 2 a (Lemma 20). 

We take a = 2s. To compute the 2e-Delaunay centers of a, we will use a pyramid data structure. 
The pyramid consists of at most log | levels. Each level h > 0 is a grid of resolution 2 ~ h z. The grid 
at level 0 consists of the single cell, D. Each node of the pyramid is associated to a cell of a grid. 
The children of a node v correspond to a subdivision into 2 d subcells (of the same size) of the cell 
associated to za The leaves are associated to the cells of the finest grid whose cells have diameter e. 

A node of the pyramid that is intersected by all the bisecting hyperplanes of <7 will be called a 
full node or, equivalently, a full cell. By our definition of W, a cell of the finest grid contains an 
element of W at its centre. The fidl-leaf-points are the elements of W associated to full cells at the 


finest level. By Lemma 22 the full-leaf-points are 2e-centers for a. In order to identify the full-leaf 
cells, we traverse the full nodes of the pyramid starting from the root. Note that to decide if a cell 
is full, we only have to decide if two corners of a cell are on opposite sides of a bisecting hyperplane, 
which reduces to evaluating a polynomial of degree 2 in the input variables. 

We will compute a bound n CT (e) on the number of full cells in the pyramid. Consider a level in 
the pyramid and the associated grid G u whose cells have diameter oj. Any point in a full cell of G u 
is at distance u from the bisecting hyperplanes of a and so, by Lemma [22] is a 2ca-center for a. By 
construction of D, the distance between a point w E Ll and a vertex p E cr is bounded by 


\w 


— p\\ < (3A 7 + 3 a)Vd < 4A 'Vd. 


Therefore, it follows from Lemma [21] that all the full cells of G w are contained in a ball of radius 
centered at c a . Letting C be a hypercube with the dimensions of a cell in G u , a volume 
argument gives a bound on n u , the number of full cells in G u : 

vol(.B) 

n-u: A 777ft = Ud 

vol(G) 
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where U d is the volume of the Euclidean ball of unit radius. 

Observe that does not depend on u;. Hence we can apply this bound to all levels of the 
pyramid, from the root level to the leaves level. Since there are at most log 5A £ V ^ levels, we conclude 
with the following lemma: 


Lemma 23 The number of full cells is less than 

. . U d {M) d hVd\' 

n ’ {£) =Wowr g —'■ 

n a (s) is also a bound on the time to compute the full cells. 


5.3 Construction of Del^L', W) 


By Lemma 20 all the simplices incident to a vertex p' of Del i£ (L 7 ) are contained in N(p') = 


L' n B(p', 2X' + 4e), and it follows from the fact that L' is a (A', //)-net that 

2 o(d) 

lA'WI < 

In the first step of the algorithm, we compute, for each p' E L', the set N(p'), and the set of 
ri-simplices 

Cd(p’) = {cr = {p'} U a : \a\ = d and d C N(p') \ {p 7 }} . 

Observe that 

(\N(p')\\ 2°^ 


' c ^' = v i j~vr 

We then extract from C d (jp' ) a subset WC d (p') of simplices that have a full-leaf-point that is a 
2e-Delaunay center, and have a number of full cells less than or equal to 


no 00 = f U d 


f 4 d \ d , 5\ZdA' 

(e^j log ^ 


This is done by applying the algorithm of Section 5.2 with a twist. As soon as a d-simplex appears 
to have more than no(e) full cells, we stop considering that simplex. The union of the sets WC d (p') 
for all p' E L 7 is a subcomplex of Del 2e (L 7 ) called Delg £ (L',W). It contains every d-simplex cr 
in Del 2e (L 7 ,lT) that has a 2e-Delaunay center in IT at a distance less than e from its actual 
circumcenter, and satisfies the thickness criterion ©(cr) > ©q. Note, however, that we do not 
claim that every simplex that has at most no(e) full cells is ©o-thick. Algorithm [3] describes the 
construction of DelQ £ (L, W). As noted above, \N(p')\ = for any p' E L\ and all the N(p') can 

be computed in 0(|L 7 | 2 ) time by a brute fore method. But assuming we have access to “universal 
hash functions” then we can use “grid method” described in | IHP111 Chapter 1] with the sparsity 
condition of L to get the complexity down to 

2°W\L'\ 

Using the facts that A 7 < 2A and ft' > see Lemma 
the algorithm is 

\L\ 
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we conclude that the total complexity of 


O 


©d d lo § £ 


and is therefore sublinear in \W\. The constant in O depends only on d. 
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Algorithm 3 Construction of Deld £ (L 7 , W ) 

input: a (A', ju')-net L', a grid W of resolution e, and ©o (defined in Eq. [ 5 ] and [6]) ) 

S := 0 

for each p’ £ L' do 

compute the sets N(p') = L' Cl B(p', 2A 7 + 4 s) and Cd{p') 

end for 

for each p' £ L' do 

for each d-simplex a £ C ( \ (p') do 

// find the simplices in Cd{p') that are in WCd{p') 

if a has a full-leaf-point that is a 2e-Delaunay center and n a (e) < no(e) then 
add a to E 

end if 
end for 
end for 

output: E the set of d-simplices of Delp £ (L 7 , W) 


5.4 Correctness of the algorithm 

We will need the following lemma which is an analog of Lemma [6] The lemma also bounds Bo¬ 
lts proof follows directly from Lemma [7j and the observation that any simplex with a protected 
circumcentre is a Delaunay simplex. 


Lemma 24 Suppose that the d-simplices in Del(L 7 ) are d-protected at their circumcenters, with 
5 = 5A 1 . If 0 < ©o < , then Del(L') C Delg £ (L / ,TT) and if, in addition, every d-simplex of 

Delp e (L',W) has a protected circumcenter, then Delo £ (L 7 , W) = Del(L 7 ). 

We first show that Algorithm [2] terminates if we deactivate the call to procedure checkQ. As 
discussed after Lemma [20] the analysis of Section [4 . 3| implies that the perturbations of Algorithm [2] 
can be expected to produce a point set L' for which all the d-simplices in Del 2e (L / , T d ) have a 

shows that we can expect the algorithm to terminate with the condition that Delg £ (L / , W) has good 
links. 

We now examine procedure check () and show that it does not affect the termination guarantee. 
By Lemma [7j if Del(L 7 ) is d-protected, then any a £ Del(L 7 ) satisfies 


d-protected circumcentre. Since this complex includes both Del(L 7 ) and Deify (L',W), Lemma 


24 


0(a) > 


8 d — 


d/T def 

24d “ 


©0- 


(5) 


Consider now a £ Delo(L 7 , W). Since the full leaves of the pyramid data-structure for a are 
composed entirely of 2e-centres at a distance less than 4\/dA 7 from any vertex of er, Lemma 
implies that, if a is ©o-thick, then 


21 


x - c n 


< 


8 \fde 
©o h'' 


This means that we can restrict our definition of Delg £ (L 7 , W) to include only simplices for which 
the set of full leaves has diameter less than • Further, we observe that if a is d-protected 
at its circumcentre, then it will have a (d — 2e)-protected full-leaf-point; this follows from the 
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triangle inequality. The check () procedure described in Algorithm [4] ensures that all the simplices 
in Delg e (Z/, W) have these two properties. It follows from the discussion above that activating 
procedure check () does not affect the termination guarantee. 

Algorithm 4 procedure check(p') 

if all d-simplices a E star(p'; Delo e (L', W)) satisfy 

1. The diameter of the full leaves is at most ■ 

2. There is a (5 — 2e)-protected full-leaf-point then 
check(p') = TRUE 

else 

check(p') = FALSE 

end if 


The fact that the algorithm terminates yields DelQ £ (L / ,lU) with good links. In order to apply 
to guarantee that Delg £ (L / ,l / L) = Del(Z/), we need to guarantee that the simplices of 


24 


Lemma 

DelQ £ (L / ,l / U) are protected. The analysis of the Moser-Tardos algorithm in Section 4.3 provides a 
possible 5-protection with 


5 = Jp, 


( 6 ) 


where J is defined in Eq. ([3]) . However, we cannot guarantee that this protection is actually obtained. 
The following lemma shows that if the perturbation p (and consequently 5) is sufficiently large, then 
the output point set will be 5*-protected for a positive 5*. 


Lemma 25 


The d-simplices in DelQ £ (L / ,lU) produced by Algorithm^ are d*-protected, with 


5* = 5- 


f 34 Vd\ 


where 5 = Jp. 


Proof of Lemma 


25 


Let a E Deln^ 


\L !, W) be a ri-simp lex produced by Algorithm]^ and let w be 
the full-leaf-point for a found by the check () procedure. Since c a is located in a full leaf-node, the 
restriction on the diameter of the leaf nodes implies 


C CT — W < 


16 \fde 
Q 0 p' 


Since w is a (5 — 2e) protected 2e-centre, we have that for all p E a and q E L' \a: 



> 

> 


> 


| q — tc|| — \\w — c 0 


\\p - k;|| + (5 - 2 e) - 


\\p ~ Co-11 + (5 - 2s) - 


16 yfds 
O 0 p' 
32 Vde 


©o P' 


\\P~Ca || + U- 


34 y/de 
©o P' 


□ 
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In order to have <5* > 0, we need a lower bound on 5, and hence on the minimal perturbation 
radius through 6 = Jp. Therefore we require: 

Jp , 34 \fde 

T^ S> l^r 

(compare with (|4)). Writing 5 = jj and ©o = , and using X' > A, we obtain the conditions under 

which Algorithm 2] is guaranteed to produce a <5*-protected Delaunay triangulation: 


Jp - 2448<P e 

— > 5 > —=-. 

4 - 5p 2 A 


The right-hand inequality is satisfied provided 


6 > 


50 d§ 
T 


We have proved 


Theorem 26 If p < p /4 and p = fl (with J defined in Eq. (HU Algorithm 2 terminates 


and outputs the Delaunay triangidation of L'. The Delaunay d-simplices are 6* -protected, as defined 
in Lemma 25, and consequently satisfy a thickness bound of 

e ( a )> s * iii/3 + n d ^ f ©,. 


The complexity of the algorithm is 


O 


8 d 


■ L l l og A 

Qd^+d 


The constants in P and O depend only on d. 


6 Boundaries and Phantom points 

We work on because our algorithm is not designed to compute a triangulation with boundaries. 
If we are given a finite set P C we scale it into the unit box, and then extend it periodically. In 
effect, by working on the torus, we are adding “phantom points” to the given point set, in order to 
avoid the boundary problem. 

The Delaunay triangulation that we compute will coincide with the Delaunay triangulation of 
the convex hull of P C only for points sufficiently far from the boundary of conv(P). Specifically, 
suppose P is a q-sampl^jof conv(P), and let L = P U Q where Q is a set of phantom points (e.g., 
periodic copies of P ) such that Q n conv(P) = 0. Let Del|(P) be the Delaunay triangulation of P 
restricted to 

D- f (P) = |x € conv(P) | d(x, <9(conv(P))) > 7 } , 

i.e., Del|(P) is the set of simplices in Del(P) that have a Delaunay ball centered in D- / (P). Then 
Del|(P) C Del(L). Indeed, by construction, a Delaunay ball B for a E Del|(P) must be contained 
in conv(P), and since Q n conv(P) = 0, it follows that B is also a Delaunay ball in L. 

1 If D C R , we say that P is a 7-sample for D if any point in D is a distance less than 7 from a point in P. 
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The extended set L C will be a A-sample for the unit hypercube, and we observe that this 
sampling radius A may be larger than 7 in general. But that is not relevant to the above observation. 
The quantity 7 is only relevant for describing Del|(P), but this complex is not considered by the 
algorithm: the algorithm is dependent only on the parameters A and p that describe L. 

Our demand that A < 1/4 is the price we pay for the convenience of working on the flat torus, 
where we identify the phantom points Q as copies of points of P. This allows us to describe the 
algorithm by considering only the given set of points in the unit hypercube. 

We now briefly describe an alternative approach which involves remaining in JR rf and explicitly 
adding additional points to the given point set P so that the boundary of the whole point set does 
not interfere with our guarantees on the simplices whose vertices are in P. 

If we are interested in using witnesses to compute the Delaunay complex on a finite vertex set 
in M d , we need to confine our considerations to a bounded domain D C M rf . But this is not possible 
if a Delaunay simplex can have an arbitrarily large circumradius. 

Assuming only that the finite set of distinct points P contains a non-degenerate d-simplex, we 
can introduce a set Q of phantom vertices such that the set P U Q is a (A, p)-net for some ball 
D C with finite radius. Without loss of generality (scaling if necessary), we can assume that 
P is contained in the unit ball centred at the origin. Let Q be the set of 2d points defined by the 
intersection of the coordinate axes with the boundary of an axis-aligned cube centred at the origin, 
and whose sides are of length 4 \fd. 

The convex hull of Q is the cross-polytope, the dual of the cube. Each facet is an equilateral 
simplex with circumradius 2\J(d — 1). Any point of P is at least unit distance from the affine hull 
of any boundary facet, and an elementary calculation shows that the circumradius of any simplex 
in Del(P U Q ) must be less than 2d — It follows then that regardless of the dimension d, the 
points of P U Q, as well as the circumcentres of all the simplices in Del(P U Q ), are contained in a 
ball D of radius 2d centred at the origin. 

Any Delaunay ball for a facet of conv(Q) that is centred at a distance greater than 2d —1/2 from 
the origin is protected in Del(P U Q), and being equilateral, these facets have maximal thickness. 
These observations imply that the guarantees of our algorithms can be shown to apply if we take 
L = P U Q, but only perform perturbations and link tests on the points in P. For the analysis in 
this case, the lemmas which reference star 2 (p) would need to be modified. Specifically, by using 
Proposition [27] directly, instead of the simplified Lemma [5j we can still assume the protection of 
lower-dimensional simplices as is required for the Witness algorithm of Section [4j Similarly, the 
thickness of Delaunay simplices that have some vertices in Q can be ensured by using the inherent 
thickness and protection of the simplices that have all vertices in Q as a substitute for the protection 
of neighbouring Delaunay d-simplices required by Lemma [7] 

Therefore, in principle the algorithm can be adapted to work with a finite set, rather than 
a periodic one, but there is no clear advantage. The algorithm, and its analysis would become 
more complicated. Also, the sampling radius A will be very large in this setting, implying that the 
sparsity parameter p would be very small, adversely affecting the running time. The subcomplex 
of Del(P U Q ) that is guaranteed to belong to Del(P), is exactly the complex Del|(P) mentioned 
above. 

If P C is a 7 -sample for conv(L > ), and we could identify the points P C P that are at a 
distance at least 4y from <9conv(P), then we could perform the algorithm by only applying the 
link test on the points of P (the deep interior points |BDG13j ). In this case, the algorithm will 
terminate, but only the simplices incident to points of P would be guaranteed to be Delaunay, and 
this is generally a smaller subcomplex than Del|(P). 
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7 Discussion 


In earlier work (BDG14| we described a perturbation algorithm to produce a d-protected point set for 
Delaunay triangulations. By employing the Moser-Tardos algorithmic framework, as we have done 
here, we obtain a much simpler analysis for the same result. Indeed, it is precisely the construction 
of a d-protected point set that allows us, in Section [I] to produce a witness complex that coincides 
with the Delaunay triangulation. In the previous work, each point is only perturbed once, whereas 
here points are only perturbed as needed, but the Moser-Tardos result ensures that only a linear 
number of perturbations are made in total. As indicated in ([3]), the new analysis yields protection 
(and therefore thickness) of the order ( l/2)°^ d2 \ improving on the (1/2)°^'^ obtained in [ BDG14] , 
The techniques described in this paper can also be applied to construct weighted, Delaunay 
triangulations and Delaunay triangulations on manifolds as will be reported elsewhere. 
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A Inheritance of protection 

If L C is d-protected, then by definition every d-simplex a G Del(L) has a d-protected Delaunay 
ball. In this section we quantify the amount of protection this guarantees for simplices of lower 
dimension. 

We will demonstrate: 


Proposition 27 Suppose L C is a locally finite set, and a G Del(L) is a non-degenerate k- 
simplex whose Delaunay balls all have radius less than A. //Vor(cr) is bounded, and every d-simplex 
in Del(L) that has cr as a face has circumradius at least \fi\, and is 5-protected, with 6 = 5 A, and 
5 < fi < 2, then a is 


(fi + 5)5 
4 (d — k + 1) 


-protected. 


Specifically, letting j = d — k, there exists j + 1 d-simplices {do, ..., aj} C Del(L), each having a 
as a face, and a is -protected at 


c = 


j + lj 


£■ 


i =0 


If k = 0, then a = {p}, and it is easy to see directly that a must be d-protected: if q G L and q p, 
then there is a d-simplex a d G Del(L) such that p G cr d and q 0 cr d . Since a d is d-protected, we 
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have ||p — g|| > 6, and it follows that the trivial Delaunay ball for p is <5-protected. Therefore the 
smallest lower bound on the protection of lower dimensional simplices happens when k = 1. With 
this observation, Lemma [5] follows immediately from Proposition |27j and the fact that the protection 
of the d-simplices in star(p) implies there are no degenerate simplices |BDG13l Lemma 3.5]. 

A standard exercise establishes that L is generic if and only if the Voronoi cell associated to any 
L’-simplex has dimension j = d — k. Rather than requiring L to be generic, Proposition |27| makes 
the more local assumption that a is non-degenerate. In this setting, this is sufficient to ensure that 
Vor(er) has dimension d — k. To see this, first observe that a short argument reveals that since all 
the d-simplices that contain cr are d-protected, none of them can be degenerate. Then another short 
argument shows that all the faces of a non-degenerate protected d-simplex must have Voronoi cells 
of the correct dimension. 

It will be convenient to work with squared distances and the “lifting map”, defined below, and 
so we introduce a modified notion of protection, called power-protection, to accommodate this. 

A.l Power-protection, and its relation to protection 

We say that a Delaunay ball B = B(c, r) for o is 5 2 -power-protected if 

||c — g|| 2 — r 2 > <5 2 , V q E L\cr. 

The following result is a special case of a result on power-protection proved in |BDGQ141 . Lemma 8]. 


Lemma 28 Suppose L C is a locally finite set, and cr E Del(L) is a non-degenerate k-simplex. 
//Vor(fj) is bounded, and every d-simplex in Del(L) that has a as a face is 5 2 -power-protected, then 
a is 

5 2 

—— — — ^ -power-protected. 

Specifically, letting j = d — k, there exists j + 1 d-simplices {cro, ..., aj} C Del(L), each having 
a as a face, and a is -power-protected at 


c = 


3 + If 


E' 


i =0 


Protection implies power-protection. Specifically, Let B = B(c,r ) be a Delaunay ball for a, and 
let R = rnirigg/-^ ||c — q\\. If B is 5- protected, then R — r > 5. Assume that the circumradius of a 
is at least A. Then 

r > 2 ^’ an d R > r + 5 > -p\ + 6. 

Thus, using 5 = <5A, we find 

R 2 — r 2 = (R + r)(R — r) > (pX + 5)(R — r) > pXS + 5 2 = ( p5 + 5 2 ) A 2 , 


which implies: 


Lemma 29 Suppose a E Del(L)) has circumradius at least pX. If a is 5-protected, then it is 
(( p5 + 5 2 )X 2 )-power-protected. 
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Going the other way, assume that a is J 2 -power-protected, so that R? —r 2 = 5 2 . Assume further, 
that r < A, and that 5 2 < 8A 2 . Then a is Y-protected, where 


S' = (R-r) = 


(R + r) 


(R 2 - r 2 ) > 


2A +Y 


If Y < 2A, then 


^4A’ 


(7) 


and by assumption we have 2A > anyway, so Inequality ([7]) is always true. We have shown: 


Lemma 30 Suppose B = B(c,r) is a 6 2 -power-protected Delaunay ball for a E Del(L). If r < A, 
and 5 2 < 8A 2 , then B is a 5'-protected Delaunay ball for a, where 

5 = 4A' 


As the notation suggests, we can employ both Lemma 29 and Lemma 30 whenever L is a (A, /i)-net. 


Proof of Proposition 21 Given the Yprotection and lower bound on the radius of the d-simplices 


28 


that are co-faces of a, Lemma 29 ensures that these simplices are ((/r + <5)<5A)-power-protected. Then 
Lemma 


ensures that a is -power-protected at c* 


We observe that the condition 5 < p, < 2, implies that (p6 + h 2 )A 2 < 8A 2 , and so the bound on 

: radii 
i+<5)(S 
4(d—k-\-l) 


the radius of the Delaunay ball for a centred at c* allows us to employ Lemma 30 to find that a is 
njrzrTTx -protected at c*. 


□ 


References 

[AEM07] D. Attali, H. Edelsbrunner, and Y. Mileyko. Weak witnesses for Delaunay triangulations 
of submanifolds. In Proc. ACM Sympos. Solid and Physical Modeling , pages 143-150, 
2007. 

[AS08] N. Alon and J. H. Spencer. The Probabilistic Method. Wiley-Interscience, New York, 
USA, 2nd edition, 2008. 

[BDG13] J.-D. Boissonnat, R. Dyer, and A. Ghosh. The Stability of Delaunay Triangulations. 

Int. J. on Comp. Geom. (IJCGA), 23(4 Sz 5):303-333, 2013. 

[BDG14] J.-D. Boissonnat, R. Dyer, and A. Ghosh. Delaunay stability via perturbations. Int. J. 

Comput. Geometry Appl ., 24(2) :125—152, 2014. 

[BDG014] J.-D. Boissonnat, R. Dyer, A. Ghosh, and S. Y. Oudot. Only distances are required 
to reconstruct submanifolds. CoRR , abs/1410.7012, 2014. Submitted to Computational 
Geometry: Theory and Applications. 

[BM14] J.-D. Boissonnat and C. Maria. The Simplex Tree: An Efficient Data Structure for 

General Simplicial Complexes. Algorithmica , 70(3):406-427, 2014. 

[Del34] B. Delaunay. Sur la sphere vide. Izv. Akad. Nauk SSSR, Otdelenie Matematicheskii i 
Estestvennyka Nauk, 7:793-800, 1934. 


21 












[dS08] 

[dSC04] 

[FKMS05] 

[HallOl 

[HP11] 
[MS 10] 

[MT10] 

[SSG89] 


V. de Silva. A weak characterisation of the Delaunay triangulation. Geometriae Dedicata, 
135(l):39-64, 2008. 

V. de Silva and G. Carlsson. Topological estimation using witness complexes. In Proc. 
Sympos. Point-Based Graphics, pages 157-166, 2004. 

S. Funke, C. Klein, Iv. Mehlhorn, and S. Schmitt. Controlled perturbation for Delaunay 
triangulations. In Proc. 16th ACM-SIAM Symposium on Discrete Algorithms (SODA), 
pages 1047-1056, 2005. 

Dan Halperin. Controlled perturbation for certified geometric computing with fixed- 
precision arithmetic. In Komei Fukuda, Jorisvander Hoeven, Michael Joswig, and 
Nobuki Takayama, editors, Mathematical Software ICMS 2010, volume 6327 of Lecture 
Notes in Computer Science, pages 92-95. Springer Berlin Heidelberg, 2010. 

S. Har-Peled. Geometric Approximation Algorithms. American Mathematical Society, 

2011. 

D. L. Millman and J. Snoeyink. Computing Planar Voronoi Diagrams in Double Preci¬ 
sion: A Further Example of Degree-driven Algorithm Design. In Proc. 26th ACM Symp. 
on Computational Geometry, pages 386-392, 2010. 

R. A. Moser and G. Tardos. A constructive proof of the generalized Lovasz Local Lemma. 
Journal of the ACM, 57(2), 2010. 

D. Salesin, J. Stolfi, and L. J. Guibas. Epsilon Geometry: Building Robust Algorithms 
from Imprecise Computations. In Proc. 5th Ann. Symp. on Computational Geometry, 
pages 208-217, 1989. 


22 



